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Nonequilibrium behavior and dynamic phase transition properties of a kinetic Ising model under 
fvj ' the influence of periodically oscillating random-fields have been analyzed within the framework of 

effective field theory (EFT) based on a decoupling approximation (DA). Dynamic equation of mo- 
C^ • tion has been solved for a simple cubic lattice {q — 6) by utilizing a Glauber type stochastic process. 

^Sj ' Amplitude of the sinusoidally oscillating magnetic field is randomly distributed on the lattice sites 

according to bimodal and trimodal distribution functions. For a bimodal type of amplitude distri- 
ct ■ bution, it is found in the high frequency regime that the dynamic phase diagrams of the system in 
temperature versus field amplitude plane resemble the corresponding phase diagrams of pure kinetic 
Ising model. Our numerical results indicate that for a bimodal distribution, both in the low and 
^h . high frequency regimes, the dynamic phase diagrams always exhibit a coexistence region in which 
fNl ' the stationary state (ferro or para) of the system is completely dependent on the initial conditions 
whereas for a trimodal distribution, coexistence region disappears depending on the values of system 
. parameters. 
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Ising model in a quenched random magnetic field (RFIM) has attracted a considerable interest in the last three 
decades. The model which is actually based on the local fields acting on the lattice sites which are taken to be random 
according to a given probability distribution was introduced for the first time by Larkin [1| for superconductors and 
later generalized by Imry and Ma ^ . A lower critical dimension d^ of the RFIM has remained an unsolved mystery 
for many years and now it is well established that a transition should exist in three and higher dimensions for finite 
temperature and randomness, which means that dc — 2 (2t6|- A great many experimental works have paid attention to 
the equilibrium properties of RFIM and quite noteworthy results have been obtained. For instance, it has been shown 
that diluted anti-ferromagnets such as FexZni_xF2 0,131, Rb2CoxMgi_xF4 I9l,ll0|, and CoxZni_xF2 [101 ^^ ^ uniform 
magnetic field just correspond to a ferromagnet in a random uniaxial magnetic field [ll|, [l2|- From the theoretical 
point of view, equilibrium RFIM has been studied by a wide variety of techniques such as mean-field theory (MFT) 
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[I3|-ll6|, effective-field theory (EFT) [17H2l|, Monte Carlo (MC) simulations [22-|25l, and the series expansion (SE) 
method [2^]. Based on these theoretical works, it is well known that different random-field distributions may lead to 
different phase diagrams and the presence of quenched randomness constitutes an important role in material science, 
since it may induce some important macroscopic effects on the thermal and magnetic properties of real materials. 

When a ferromagnetic material is subject to a periodically varying time dependent magnetic field (kinetic Ising 
model) , the system may not respond to the external magnetic field instantaneously which causes interesting behaviors 
due to the competing time scales of the relaxation behavior of the system and periodic external magnetic field. At high 
temperatures and for the high amplitudes of the periodic magnetic field, the system is able to follow the external field 
with some delay while this is not the case for low temperatures and small magnetic field amplitudes. This spontaneous 
symmetry breaking indicates the presence of a dynamic phase transition (DPT) [27| | which shows itself in the dynamic 
order parameter (DOP) which is defined as the time average of the magnetization over a full period of the oscillating 
field. DPT properties of kinetic Ising model have been firstly observed theoretically by Tome and Oliveira within the 
framework of MET [28|] . Since then, much attention has been devoted to investigate the dynamic nature of the phase 
transitions by means of several theoretical and experimental works. On the theoretical side, non-equilibrium phase 
transition properties of kinetic Ising model have been widely investigated by various techniques [29|-|36|. Besides, on 
the experimental picture, a DPT occurring for the high frequency magnetic fields was studied by Jiang et al. [37| using 
the surface magneto-optical Kerr-effect technique for epitaxially grown ultrathin Co films on a Cu (001) surface. Eor 
a [Co(4A°)/Pt(7A°)] multi-layer system with strong perpendicular anisotropy, an example of DPT has been observed 
by Robb et al. |38| . They found that the experimental non-equilibrium phase diagrams strongly resemble the dynamic 
behavior predicted from theoretical calculations of a kinetic Ising model. It is clear from these works, there exists a 
strong evidence of qualitative consistency between theoretical and experimental studies. 

On the other hand, non-equilibrium stationary states and dynamic phase transition properties of REIM have not 
been well understood yet, and there exists a limited number of works in the literature [39|-|43|. Eor instance, Hausmann 
et al. [39| considered the behavior of an Ising ferromagnet under the influence of a fast switching, random external 
field. According to the analytical results based on MET for the stationary state of the system, they observed a novel 
type of first order phase transition which has also been verified by their extensive MC simulations. Paula et al. [40] 
determined the stationary states of REIM by using MET and constructed the phase diagrams from the stationary 
states of the magnetization as a function of temperature and field amplitude. They found that the continuous phase 
transitions coincide with the equilibrium ones 'l5|, while the first-order transitions occur at fields larger than the 
corresponding values at equilibrium. In addition, they also observed that the difference between the fields at the 
limit of stability of the ordered phase and that of the equilibrium is maximum at zero temperature and vanishes 
at the tricritical point. Eurthermore, Acharyya [4l| studied the non-equilibrium dynamic phase transition, in the 
two-dimensional kinetic Ising model in the presence of a randomly varying magnetic field both by MET and MC 
simulations and discerned that in contrast to the results found in Ref. |39|, the nature of the transition is always 
continuous. In a recent work, Crokidakis [42j performed MC simulations on cubic lattices for a non-equilibrium Ising 
model that stochastically evolves under the simultaneous operation of several spin-flip mechanisms where the local 
magnetic fields change sign randomly with time due to competing kinetics. Erom the numerical results, it has been 
predicted that there exist first-order transitions at low temperatures and large disorder strengths, which correspond 
to the existence of a non-equilibrium tricritical point at finite temperature. Very recently, Costabile et al. [43| have 
studied the dynamical phase transitions of the kinetic Ising model in the presence of a random magnetic field by 
using EET with correlations where the EET dynamic equation has been given for the simple cubic lattice {q — 6), 
and the dynamic order parameter has been calculated. It has been observed that the system presents ferromagnetic 
and paramagnetic states for low and high temperatures, respectively. Apart from this, they have predicted a non- 
equilibrium tricritical point in a phase diagram in the temperature versus applied field amplitude plane. They have 
also compared the results with the equilibrium phase diagram [17lll8| . where only the first-order line is different. In the 
theoretical works mentioned above, the random field effects have been taken into account either by a given probability 
distribution function (random in space) namely a bimodal distribution, or by generating a new configuration of random 
fields uniformly at each time step (random in time) . 

In the present paper, we have studied dynamic phase transitions and stationary states of REIM driven by a 
periodically varying time dependent magnetic field on a simple cubic lattice. Amplitude of the applied magnetic field 
is sampled from both bimodal and trimodal probability distributions. Eor this purpose, we organized the paper as 
follows: In Sec. |lT]we briefiy present the formulations. The results and discussions are summarized in Sec. IIIIl and 
finally Sec. IIVI contains our conclusions. 



II. FORMULATION 

We consider a three dimensional Ising ferromagnet ( J > 0) defined on a simple cubic lattice with a time dependent 
external magnetic field. The time dependent Hamiltonian describing our model is 



n=-jY. S,Sj -Y,H^{t)S^. 



(1) 



<ij> 



where the first term is a summation over the nearest-neighbor spins with S'i = ±1 and Hi{t) is a time-dependent 
external oscillating magnetic field which is given by 



Hi{t) = Ho^cos{ujt), 



(2) 



where Hoi is the amplitude of the external magnetic field acting on the site i, and w denotes the angular frequency 
of the oscillating external field. The amplitude of the field is distributed according to a given probability distribution 
function. The present study deals with a trimodal field distribution which has a form 



P{Ho^) = p5{Hoi) 



l-p 



[6{Ho, - Ho) + 8{Hn + Ho)] 



(3) 



In order to describe the dynamical evolution of the system, we follow a Glauber type stochastic process ^] . The 
dynamical equation of motion can be obtained by using the master equation as follows 



dt 



-(S',)-h/tanh 



E, + HXt) 
keT 



(4) 



where t is the transition rate per unit time, Ei — J^j Sj is the local field acting on the lattice site i, and ks and T 
denote the Boltzmann constant and temperature, respectively. 

In Eq. Q , we set r at unity. If we apply the differential operator technique [4^, |40| in Eq. ^ by taking into 
account the random configurational averages we get 



dm 
'dt 



/ q—Q 



-m - 



Yl cosh( JV) + S^ sinh( JV) ) ) F(a;)U=o, 



(5) 



where m = {{Si))r represents the average magnetization, V — d/dx is a differential operator, q is the coordination 
number of the lattice, and the inner (...) and the outer (...)r brackets represent the thermal and configurational 
averages, respectively. The function F{x) in Eq. ([5|) is then defined by 



F{x) 



dHo,PiHo,) tanh 



x + H,(t) 



knT 



(6) 



When the right-hand side of Eq. ([S]) is expanded, the multispin correlation functions appear. The simplest approxi- 
mation, and one of the most frequently adopted is to decouple these correlations according to 



{{S.S,...Si))^^{{S,))^{{S,))^...{{Si))^, 



(7) 



for i 7^ j ^ ... 7^ / [47]. If we expand the right-hand side of Eq. ([5]) within the help of Eq. ([7]) then we obtain the 
dynamical equation of motion as follows 



dm 
~dt 



q=6 



(8) 



The coefficients in Eq. (|5]) are defined as 



A, 



29 



EE 



r 



{-!)'' exp [{q ~2r- 2s) JV] F{x) 



x=0, 



j = 0,l,...,g. 



(9) 



These coefficients can be calculated by employing the mathematical relation exp{aV)F{x) — F{x -I- a). Eq. ([5]) can 
be regarded as a kind of initial value problem and the solution can be easily found benefiting from the initial value of 



the average order parameter mo by using the fourth order Runge-Kutta method (RK-4) . For selected values of the 
Hamiltonian parameters and temperature, the time dependence of magnetization converges to a finite value after some 
iterations i.e. the solutions have property m{t) = m{t + 2tt/uj) for arbitrary initial value of the magnetization (toq). 
Thus, by obtaining this convergent region after some transient steps (which depends on Hamiltonian parameters and 
the temperature) the DOP which is the time average of the magnetization over a full cycle of the oscillating magnetic 
field can be calculated from 

Q = ^im{t)dt, (10) 

where m(t) is a stable and periodical function which can be one of the two types, according to whether it has the 
following property or not [28| 



m{t) = —m{t + tt/cj). 



(11) 



A solution satisfying Eq. ()11|) is called symmetric solution which corresponds to a paramagnetic (P) phase where 
the magnetization oscillates around zero whereas the solution which does not ensure Eq. (jlip is called nonsymmetric 
solution, and it corresponds to ferromagnetic (F) phase where the magnetization oscillates around a non-zero value. 
In these two cases, the observed behavior of the magnetization is regardless of the choice of the initial value of 
magnetization. On the other hand, in contrast to the equilibrium RFIM, there exist coexistence regions (F+P phases) 
in the phase diagrams in temperature versus field amplitude plane where the stationary state of non-equilibrium RFIM 
problem depends on the initial value mo of the time dependent magnetization. Furthermore, it is not possible to obtain 
the free energy for kinetic models in the presence of time dependent external fields. Hence, in order to determine the 
type of the dynamic phase transition (first or second order) , it is convenient to check the temperature dependence 
of DOP. Namely, if the DOP decreases continuously to zero in the vicinity of critical temperature, this transition is 
classified as of the second order whereas if it vanishes discontinuously then the transition is assumed to be of the first 
order. 

III. RESULTS AND DISCUSSION 

In this section, we discuss how the random-fields effect the phase diagrams of the kinetic Ising model. Also, in 
order to clarify the type of the dynamic phase transitions in the system, we give the temperature dependence of the 
dynamic order parameter. 
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FIG. 1: Dynamic phase diagram of the pure kinetic Ising model in (fesT/oJ — Ho/qJ) plane for oscillation frequency values 
oj = 0.5 (left panel) and cj — 1.0 (right panel), for comparison with Refs. [351,1391 . Solid (dashed) lines correspond to second 
(first) order phase transitions and solid symbols represent dynamic tricritcal points. 



A. Non-equilibrium phase diagrams of the kinetic model: pure case 

In order to provide a testing ground for our calculations, we have primarily studied the dynamic phase diagrams 
of the kinetic Ising model under an oscillating magnetic field where the amplitude of the externally applied field 
was taken as a uniformly constant value. This model defines non-equilibrium properties of the pure system and has 
been examined previously within the framework of EFT [3a, [Sg. In these works, the authors investigated the phase 
diagrams of the system in a [hgTc/qJ — Hq/qJ) plane where q is the coordination number of the lattice. From these 
works, we see that, for the oscillation frequency values uj = 0.5 and 1.0, the location of dynamic tricritical (DTC) 
point was identified imprecisely. Hence, in order to compare our results with those found in Refs. [35|, |36[, we depict 
the phase diagrams in Fig. [T]in the same plane. It is clear from Fig. [l]that our numerical values of DTC point for 
w = 0.5 and 1.0 agree well with [36| whereas they are quite different from those obtained in Ref. [35| which is probably 
due to the fact that in Ref. [35| , insufficient number of data points was used to construct the phase diagrams of the 
system. 

B. Non-equilibrium phase diagrams of RFIM for a bimodal distribution 

The distribution function given in Eq. ([3]) corresponds to a bimodal field distribution for p = where the amplitude 
of the oscillating field can be either +Ho or —Ho with equal probability. In this case, the system can be thought as 
a spin system under the influence of two oscillating external field sources. In Refs. [40, 143 a similar model has been 
studied within MFT and EFT, respectively where the authors did not consider any oscillating external magnetic field. 
Main conclusion of those studies was that the dynamic second order phase transition lines in (kBTc/J — Hq/J) plane 
coincide with the equilibrium counterparts [15J, llTI . Il8| whereas the maximum and minimum differences between the 
dynamic and equilibrium first order phase transition lines were observed at the zero temperature and at the tricritical 
point, respectively. However, as seen in Fig. [5^, if the external magnetic field has an oscillatory character with an 
amplitude which is applied at random then a completely different situation arises. Namely, the second order phase 
transition fields of non-equilibrium RFIM with oscillating random external fields are greater than those obtained for 
static RFIM [IJ, [lM|- On. the other hand, according to the dynamic phase diagram shown in Fig. [5^, the system 
exhibits F phase for high temperature and low field amplitude values. In addition, as the oscillation frequency uj 
increases then DTC point depresses and the F-|-P phase region gets narrower. In Fig. [2J5, variations of DTC point 
coordinates H^/J and ksT* / J with respect to the frequency w are plotted. As seen in this figure, Hq/J grows whereas 
UbT^ / J decays with increasing a; and these values saturate at Hq/J = 4.512 and ksT^ / J — 1.54. This means that 
the coexistence region in the phase diagrams gets narrower but does not disappear with increasing uj. In other words, 
for sufficiently high frequencies, both the location of DTC point and the area of F-l-P region become independent of 
oscillation frequency of the external field. We also note that F region is always independent of frequency for the whole 
range of w. 

In Fig. [31 we show the temperature dependence of the dynamic order parameter Q as a function of the oscillation 
frequency w corresponding to the phase diagrams depicted in Fig. [2l As we mentioned before, increasing or decreasing 
frequency values does not affect the area of the F region. Namely, for a fixed Hq/J value, the system always undergoes 
a dynamic phase transition on the F-P phase boundary line, and the critical temperature value is independent of to. 
This situation is independent of the selection of the initial magnetization for Hq/J < Hq/J whereas it is valid for 
Hq/J > Hq/J only if one starts with an initial condition ttiq = 0. As is clearly seen in Fig. [31 all curves coincide with 
each other for Hq/J = 2.0. The inset shows the average magnetization m{t) as a function of time t at ksT/ J = 4.0 
and for oj = 0.1,0.3 and 0.5 where we see that m(t) curve oscillates with smaller amplitude as w increases, but the 
average value, i.e. DOP (Q) over a complete cycle of the magnetic field does not change. On the other hand, as shown 
in the right panel of Fig. [31 for sufficiently strong amplitudes of the external field, such as Hq/J = 4.0, DOP versus 
temperature curves exhibit different characteristics with increasing frequency, although the critical temperature does 
not change, and interestingly, for sufficiently high frequency values such as w > 0.5 DOP curves resemble those of the 
pure kinetic Ising model driven by a periodic external field with the high oscillation frequency. 

In Fig. [31 dynamic phase diagrams of the system are depicted for a wide range of oscillation frequency values. 
In this figure, we observe that DTC point depresses for a while and the area of the coexistence region in the phase 
diagrams slightly gets narrower and then remains unchanged with increasing frequency values. The low frequency 
phase diagrams of the system are found to be completely different from those of the pure kinetic Ising model [3g| 
whereas for the high frequency regime, dynamic phase diagrams of the system are exactly identical to those of the pure 
case. These observations originate from the symmetry of the random field distribution. Namely, for a bimodal field 
distribution, we have two oscillating magnetic field sources. Initially, half of the lattice sites are under the influence 
of a periodically oscillating field (source-1) with an amplitude Hq/J, and the remaining spins on the other half of the 
lattice are influenced by an oscillating external field with an amplitude —Hq/J (source-2). As the time progresses, 
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FIG. 2: (Color online) (a) Phase diagrams of the kinetic RFIM in (kBTc/J — Ho/J) plane for a bimodal distribution with uj = 0.1 
and 0.5, in comparison with the static model [17, 18]. Solid (dashed) lines correspond to second (first) order phase transitions 
and solid symbols represent tricritical points, (b) Variation of the coordinates Hq/J and ksT^ / J of dynamic tricritical point 
as a function of frequency ui. 
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FIG. 3: (Color online) Temperature dependence of the dynamic order parameter Q as a. function of frequency u for Ho/J — 2.0 
(left panel) and 4.0 (right panel). Insets in both panels show the time variation of the average magnetization m{t) at ksT/J = 
4.0 (left) and 1.75 (right). 



due to the sinusoidal form of the periodic fields, magnetic field sources acting on the lattice sites change their sign 
instantaneously at the end of the one half of the oscillation period. Consequently, the presence of a phase difference 
between the source-1 and source-2 stimulates some unusual effects on the system. Furthermore, for sufficiently high 
frequencies of the oscillating external fields, relaxation time of the system becomes greater than the oscillation period 
of the external field. As a result of this, source-1 and source-2 fields receive identical responses from the system which 
perceive the field sources as a single oscillatory field source, hence the dynamic phase diagrams of the system resemble 
the high frequency phase diagrams of the corresponding pure model. We may also note that F+P regions in the 
ksTc/J — Hq/J phase diagrams always exist for the whole range of frequency values. 

In Fig. [51 we represent some examples of typical DOP versus temperature profiles corresponding to the phase 
diagrams shown in Fig. 2] with w = 1.0. Namely, the system always undergoes a second-order dynamic phase 
transition for Hq/J — 4.4. In this case, dynamic nature of the phase transition and the stationary state of the 
system are independent of the initial magnetization mg. On the other hand, two successive first order dynamic phase 
transitions (i.e., a first-order reentrant phenomena) occur at the value Hq/J — 4.52 with the initial condition mo = 0.0. 
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FIG. 4: (Color online) Dynamic phase diagrams of the system in a {ksTc/J — Hq/J) plane for a bimodal distribution with 
some selected values of frequency u. The solid (dashed) lines correspond to second (first) order phase transitions and solid 
symbols represent tricritical points. 

In addition, for Hq/ J = 4.6 with ttiq = 1.0, DOP curve exhibit a discontinuous jump at a phase transition temperature 
which puts forward the existence of a first-order transition on the system. For the latter two examples, dynamic nature 
of the phase transition and the stationary state of the system strictly depends on the initial magnetization. 
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FIG. 5: (Color online) Variation of the dynamic order parameter Q as a function of the temperature for a; = 1.0 with some 
selected values of Hq/J- The dashed lines correspond to first order phase transitions. 



C. Non-equilibrium phase diagrams of RFIM for a trimodal distribution 

For a trimodal field distribution of magnetic fields defined in Eq. ([3]), we plot the dynamic phase diagrams of the 
system in a {ksTc/J — Hq/J) plane in Fig. [5] for w = 0.5 and with some selected values of disorder parameter p. As 
shown in Fig. [S^, the system exhibits DTC behavior for a relatively narrow range of p. The temperature coordinate 
UbT* / J of DTC point decreases as p increases whereas the field amplitude part Hq/J increases, and consequently 
DTC point depresses to zero at a certain value of p. After the destruction of the DTC point, all dynamic phase 
transition processes are found to be of the second order. We also note that at the zero temperature, there exist a 
critical value p* = 0.53 below which the system exhibits a dynamic phase transition at a critical magnetic field. This 
value of p* is independent of w. Moreover, according to Fig. |6)d, as p increases then the area of F+P region in the 
dynamic phase diagrams gets narrower and after a specific value of p, we cannot observe any coexistence region in 
(fcsTc/ J - Ho/ J) plane. 
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FIG. 6: (Color online) Dynamic phase diagrams of the system in a (ksTc/J — Ho/ J) plane for a trimodal distribution of 
random fields with w = 0.5 and for some selected values of p. The solid (dashed) lines correspond to second (first) order phase 
transitions and symbols represent tricritical points. 
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FIG. 7: (Color online) Variation of DOP (Q) as a function of the amplitude of the oscillating field Hq/J for a trimodal 
distribution corresponding to the dynamic phase diagrams depicted in Fig. [6] with some selected values of p for mo = 0, uj — 0.5 
and for the reduced temperatures (a) ksT/J = 0.5, (b) ksT/J — 2.5. 



Variation of DOP (Q) as a function of the amplitude of the oscillating field Hq/J for a trimodal distribution 



corresponding to the dynamic phase diagrams depicted in Fig. |6] is plotted in Fig. [7] with some selected values of 
p. In this figure, initial value of the average magnetization is selected as mg = 0.0. According to Fig. [71 DOP 
curves undergo first and second order dynamic phase transitions at low and high temperatures, respectively. At 
low temperatures, critical amplitude value H^/ J at which a dynamic phase transition occurs depends on the initial 
magnetization ttiq. Namely, for ttiq — 0.0 critical Hq/ J value is greater than that obtained for toq = 1.0 which is due 
to the presence of F+P phase in the system at a fixed value of p (see Fig. |6Jd). 

Finally, for a fixed oscillation frequency a; = 0.5, let us investigate the time dependence of the average magnetization 
m{t) corresponding to the phase diagrams shown in Fig. [SJd. In Fig. |5^, we plot the time series of m(t) curves for 
p = 0, ksT/J = 1.0, and with two values of field amplitude Hq/ J. From the upper panel in Fig. |8^, we see that the 
dynamic nature of the stationary state of the system is independent of the initial magnetization toq for Hq/ J = 4.4. 
Namely, for this set of system parameters, average magnetization of the system oscillates around a nonzero value (F 
phase) after some transient time. However, for H^/ J — 4.58 (while the other parameters are the same) stationary 
state of the system (F or P) depends on the initial magnetization mo which indicates that the dynamic phase diagrams 
exhibit a coexistence region (F+P phase) in {ksTc/J — Hq/ J) plane. On the other hand, as seen in Fig. |8)d which 
is plotted for p = 0.15, we cannot observe any coexistent phase region in the system. Namely, for Hq/ J = 5.1 with 
ksT/J = 1.0, Tn{t) curves oscillate around a finite nonzero value, whereas for Hq/ J — 5.22 and ksT/J — 0.25 the 
curves oscillates around zero which means that the system is in P phase. Hence, we see that the stationary state of 
the system for a trimodal distribution of the oscillating field amplitude may be independent of ttiq, depending on the 
value of distribution parameter p. 
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FIG. 8: (Color online) Time dependence of the average magnetization m{t) for u} — 0.5 with some selected values of the 
parameters p, knT/J, and Ho/J. 



IV. CONCLUDING REMARKS 



In conclusion, we have investigated the kinetic behavior of a spin-1/2 Ising model on a simple cubic lattice {q = 6) 
under the influence of a quenched random magnetic field which oscillates periodically in time by means of an effective- 
field theory based on a standard decoupling approximation and the time evolution of the system has been presented by 
utilizing a Glauber type stochastic process. For bimodal and trimodal field distributions, we have studied the global 
dynamic phase diagrams of the system in a {ksTc/ J — Hq/ J) plane and variation of the dynamic order parameter 
with respect to the temperature and random amplitude of the oscillating field, as well as the time series of the average 
magnetization curves. We found that the dynamic behavior of the random field Ising model with a periodically 
oscillating magnetic field exhibit quite different characteristics in comparison with the static RFIM problem. 

Our numerical analysis clearly indicates that such field distributions may lead to a number of interesting phenomena. 
For example, we have found that for a bimodal distribution with sufficiently high frequency values, DOP curves and 
non-equilibrium phase diagrams of the system resemble those of the pure kinetic Ising model driven by a periodic 
external field with the high oscillation frequency. On the other hand, we have also observed that F-l-P regions in 
(ksTc/J — Hq/ J) phase diagrams always exist for the whole range of frequency values for a bimodal distribution 
of oscillatory field amplitude whereas for a trimodal distribution, area of F-l-P region gets narrower with increasing 
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distribution parameter p, and if we increase p further then we cannot observe any coexistence region in (ksTc/J — 
Hq/J) plane. 

EFT method takes the standard mean field predictions one step forward by taking into account the single site 
correlations which means that the thermal fiuctuations are partially considered within the framework of EFT. Although 
all of the observations reported in this work shows that EFT can be successfully applied to such nonequilibrium systems 
in the presence of quenched disorder, the true nature of the physical facts underlying the observations displayed in 
the system (especially the origin of the coexistence phase) may be further understood with an improved version of the 
present EFT formalism which can be achieved by attempting to consider the multi site correlations which originate 
when expanding the spin identities. We believe that this attempt could provide a treatment beyond the present 
approximation. 
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